hep-ph/0407050 



Bd — Bd mass difference in Little Higgs model 



S. RAI CHOUDHURY^*, NAVEEN GAUR^^ ASHOK GOYAL^*, Namit Mahajan^^ 

^ Department of Physics & Astrophysics, 
University of Delhi, Delhi - 110 007, India. 

^ Harish-Chandra Research Institute, 
Chattnag Road, Jhunsi, Allahabad - 211019, India. 



Abstract 

An alternate solution of hierarchy problem in the Standard Model namely, the Little Higgs model, 
has been proposed lately. In this work B'^ — B'^ mass difference in the framework of the Little Higgs 
model is evaluated. The experimental limits on the mass difference is shown to provide meaningful 
constraints on the parameter space of the model. 
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Note : In our earlier version of the paper a mixup of two symbols in our numerical computation lead 
to incorrect values of the — mass difference. A paper that appeared subsequently by Buras et.al. 
[1], presented the correct calculation of B^ — B^ mass difference in Little Higgs model. In this version 
we have corrected our numerical mistake. Our results for B^ — Bd mass difference now agrees with the 
results of Buras et.al. [1]. Some typos in our earlier version have also been corrected. 



Our understanding of the Standard Model (SM) is plagued by a major issue, the "hierarchy prob- 
lem", arising out of the enormous difference between the electroweak and the Plank scale. For quite 
some time, Supersymmetry had provided an elegant framework for solving this problem although to- 
date there is no compelling experimental evidence in its support. During the last two years, an alterna- 
tive possibility has been introduced in the literature where the Higgs mass remains small by virtue of it 
being a Goldstone boson of a global symmetry which is broken at a scale above the electroweak scale. 
These models are generically called the "Little Higgs" models and the simplest of these, the 'Littlest 
Higgs"(Ui) model [2], has the least number of additional particles involved. 

In the gauge sector, the LH model contains weakly coupled gauge bosons with masses in the TeV 
scale in addition to the SM and Z [3,4]. These mix amongst themselves causing modification of 
SM gauge couplings of W^, Z with fermions and among themselves. In the quark sector, a vector- like 
heavy top quark comes into play with mass in TeV range, which has trilinear coupling with SM gauge 
bosons. Once again the heavy top quark has mixing possibility with the SM top quark, resulting in 
modification of coupling structure of quarks with and Z. In addition, the model has charged Higgs 
bosons which introduce scalar couplings with quark. Also, a heavier photon with mass in the TeV range 
emerges, which couples both to leptons and quarks. 

The presence of these new particles as well as changes in the SM interaction vertices, can cause 
changes in a variety of measurable parameters. Some of them have already been calculated in the 
literature [3,5-8]. These results provide good constraints on the parameters entering the LH model. 
Direct experimental confirmation of several aspects of LH, e.g., the masses of the heavy t-quark and the 
doubly charged Higgs, would require sharper estimates of the parameters of the theory. It is desirable 
therefore, to work out the consequences of the LH-model for as many observable quantities as possible 
in order to sharpen the constraints on the parameter space of such a model. In this note, we report on a 
calculation of Bq — Bq and Kq — Kq mixing in the context of LH model. 

In SM, there is one basic box diagram responsible for generating the effective Hamiltonian for the 
mixing of Bq — Bq and Kq — Kq. In LH, there are many more box diagrams (as shown in Figure 1) to 
be evaluated. The couplings and propagators required for calculating these diagrams are listed in [3]. 

The effective Hamiltonian resulting for the graphs in Fig. 1 has the structure : 

Ueff = j^M^Sg{qd)v-Aiqd)v-A (1) 

with q — b,s for (Bq — Bq) and (Kq — Kq) mixing respectively. The invariant function Sq has the 
following form: 

Sq^S^^ + S^"" (2) 
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where in both Sb and Ss, the first term represents the SM contribution along with QCD corrections 
which are given in detail in [9]. The second term gives the LH contribution to the mass difference. As 
these are the corrections to the SM contribution, we do not consider QCD corrections to them which 
would arise from gluonic loops added to the diagrams of Fig 1 . The effective Lagrangian in the LH 
model to order p-^ is well approximated by : 



C 



eff - YGtI^-'^-J 



AJ 



Ml,SfQ{AJ = 2) 



where J = B, S and j = b,s for — B^ and — respectively. They are given as 



(3) 



g(Afi = 2) = {b^da)v-Aihd^)v-A 
Q{AS = 2) = iSada)v-A{Sf}dp)v-A 



(4) 



and 



qLH 
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E X^E{x„Wl,<I>)+ E KX,E{x„Xj,WL,<^) 



i=u,c,t 



(5) 



where ^i, ^ij and functions E axe defined in appendix A. Xi = mf/ rriy^, A's are the CKM factors defined 
as Aj = VidV*^ {i = u, c, t) and Ay = j^Xt and Vi/s are the CKM matrix elements. 



r\/\/\/\/x/i 



u,c,t,T 



u,c^t,T 



u,c,t,T 



u,c,t,T 



Figure 1 : Box diagrams in LH. 



We note that despite the occurrence of spinless Higgs couplings to quarks, the ultimate structure of 
the effective Hamiltonian in LH retains the same {V — A) form as in SM to order (^j^ . Given the form 

'/ is the scale at which the global SU(5) symmetry is spontaneously broken via a vacuum expectation value which is 
expected to be in the TeV range and roughly of the order of masses of heavy bosons and v is the vev of standard model 
Higgs 
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of the effective Hamiltonian, we can proceed exactly as in SM and calculate its matrix element between 
Kq — Kq or Bq — Bq states using the vacuum saturation approximation. There are no divergences in 
the SM amplitude because of the unitarity of the CKM matrix ; this statement holds even in LH modeP 
where once again the unitarity of CKM ensures that all divergences vanish to order {v/fY. Neglecting 
QCD corrections and long distance contributions we can get the mass difference to be : 

AM{Bo - Bof" = ^MBM^f^S^"" (6) 



AM{Ko - Kof"" = ^MkMUIS^" (7) 



and 

AMiKn - Kn\^^ = 

6' 
IT 

where Mb,k, fB,K are the masses and decay constants of B and K mesons respectively ^. 

It should be mentioned that the renormalization group evolution of the matrix elements has been the 
subject of much work and has been summarized in [10] and is far from trivial since the matrix elements 
are controlled by long distance dynamics and are generally parameterized by a "Bag factor" Bk,b- 
However for the neutral B meson case, the long range interactions arising from the intermediate virtual 
states are negligible because of the large B mass, being far from the region of hadronic resonances. 

The LH involves not only heavy vector bosons and quarks but also a large number of parameters 
over and above those in the SM. The global symmetry in the theory is broken at TeV range scale 
As(A5 = Anf)\ the scalar bosons, doublets and triplets, acquire vacuum expectation values v and v' 
respectively at the EW-scale, providing the convenient small parameters v/ f and v' / f . The mixing 
of the charged and neutral vector bosons results in two mixing angle parameters 9 and 6' (with c = 
CosO, s — Sind, c' — CosO' and s' — SinO'). Finally the Yukawa coupling of the fermions involves 
two parameters Ai and A2 with the combination xl = x'^+x'^ occurring frequently. To the leading order 
in {v/ f), the masses of all the heavy particles in LH can be expressed in terms of SM masses niw and 
mz as : 

mw SC V 

rrit A1A2 V 

rriH V 

The coupling of all heavy particles to SM particles as well among themselves are expressible in 
terms of these parameters with the SM ones. 

The parameter space is obviously too large. Requiring that the heavy particles have masses in TeV 
range results in the condition ^ < 10. There is another restriction arising out of the requirement that 
the mass of the triplet scalars be positive definite [3]: 



^the CKM matrix is unitary in LH upto order v"^ / p 

^The QCD corrections to these have been worked out in literature 
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We have varied f // in the range to 0.1. s, s' in range 0.2 to 0.8 and xl in range 0.2 to 0.8 in our 
numerical analysis. Other parameters used are given in the Appendix B. 

Our results for the — case are shown in Figure 2. Varying s' doesn't significantly change 
our conclusions. The corresponding — results too have similar trend. However, since there are 
large error bars in them because of QCD corrections involved, it makes it difficult to draw any definitive 
conclusions. Hence we haven't shown them. In the plots shown in Fig. 2 the shaded area corresponds 
to the mass difference AAf {B^ — Bd) = 0.5 ± 0.01 ps^^ which is consistent with current experimental 
bounds [11]. 

From Fig 2, it is easy to note that for low values of xl, there is a region of parameter space (in 
terms of the parameters s and v/ f) that is consistent with the experimental limits. Very specifically for 
Xl = 0.2, the bound on the scale / can be very close to 1 TeV for almost all the s values. However, 
as is increased, the LH contribution starts deviating significantly from the SM results. Therefore, 
in principle the experimentally allowed band for the B^ — B^ mass difference provides significant 
constraints for the parameter space of generic Little Higgs models, and in particular the Littlest Higgs 
model. 

It will be fruitful to compare the limits on the parameters coming from precision electroweak data 
[12] . To this end we recall that the Littlest Higgs model does not have the custodial symmetry inherently 
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built into it and can therefore, in principle, lead to large corrections, arising both from heavy gauge 
boson exchange diagrams and the triplet VEV. A naive way out would be to have the extra gauge boson 
masses raised by some means. However, this would spoil the motivation of circumventing the hierarchy 
problem and would also bring in the issue of fine tuning. It has been found that global fits to precision 
electroweak data imply the following bound (at 95% C.L.) on the scale / for any generic coupling 
(specifically varying c and c' between 0.1 and 0.995): 

/ > 4 TeV (9) 

It is worth noting that this very stringent bound is almost (practically) independent of any variation of 
Higgs mass upto 200 GeV. Further, the bound still holds for any order unity value of the parameters 
encapsulating the physics due to proper UV-completion of the theory. To be noted is the fact that these 
are very strong limits and the origin can simply be traced back to the absence of custodial symmetry. 
In the second reference of [12] it was noted that considering only precision electroweak data allows 
for a small region in parameter space where the bound on the scale / can be lowered to about 1 TeV. 
However, electroweak data combined with Drell-Yan production excludes this region and a combined 
analysis implies a bound very similar to the one quoted above. Constraints from low energy precision 
data like {g — 2)^ and atomic weak charge of Cesium also indicate similar bounds, though it should be 
remembered that the rather small {g — 2)^ corrections may not serve to put any meaningful bounds. 

Interestingly enough, in almost all the variations of the Littlest Higgs model [13], the constraints 
remain quite strong and generically very similar to the minimal version, though for some very specific 
choices of the parameters the constraints on / are relaxed to 1 — 2 TeV. This can be understood as 
arising due to small mixing between the two sets of gauge bosons and also small coupling between the 
fermions and heavy U{1) gauge boson. Nevertheless, these arise only in very specific models and for 
very special choices of the parameters and are not a generic feature of LH models. It may also be useful 
to keep in mind that the positivity of triplet mass squared imposes severe constraints on the triplet VEV 
and therefore the parameters or combination of parameters entering the mass squared relation. This 
strong constraint considerably reduces the allowed parameter space. Therefore, it is natural to expect 
that in the variation of the minimal model where there is no triplet Higgs, the bounds are partially 
relaxed as is the case in models which have custodial symmetry built into them. 

Turning to contribution to AM{Ba — Bo), we would indeed get a bound similar to eqn (9) above 
if we require that the LH contribution be no more than the experimental band for AM{Bd — Bd)- 
However, as has been pointed out by Buras et.al. [1] that there is a hadronic uncertainity of about 10% 
in calculation of AM{Bd — Bd). In view of this a reasonable constraint on {v/ f) could be given if the 
variation in the mass difference could be more than 10%. For a contribution of about 10% , the LH 
model would require v/ f < 0.2, which is not very useful in view of the stronger constraint like eqn(9) 
above. However, should it become possible for hadronic uncertainities to be reduced by a factor of 2 
or 3, then the bound on (f//) values becomes much lower, leading to constraints on the value of / 
comparable with the value in equation (9) above. 
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A Loop Functions 

Wl, Wh in eqn (5) refer to the light & heavy W-boson in LH; c,s are the mixing angles in LH. Various 
^'s are 

= 2c'{c'-s')Xl i^u,c (A.1) 
= 2{c^(c^-.^)+4}A? (A.2) 



Cit = 2|c2(c2-s2) + ^|a,A, (A.3) 

CiT = -xlXiXt (A.4) 

Ctr = -xiA? (A.6) 

e1 = -if/vrX^ (A.7) 

^ij = -if/vf\, (A.8) 



the functions (E) used in eqn. (5) are : 

7—1 / TT7-\ "^i'^yoIX 3 X 3 X li/X 



Xj -Xi [A 2(1- Xj) 4(1- Xj^ 
3^^_XiXj__ 

^A{l-x,){l-xj) ^ ■ ^ 

= 4(^)^^^^(^^)-^'{^i(r^-^(r^ '^-''^ 

7 rr^ — 77 U - I 1 + ——] + — ^ \ logixA 

XWh [Xj - Xi){l - Xj){l - ^) [ V J J 

3 1 

^ / , X / ^ / ^-logixwn) (A.11) 
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E{xi,WL,WH) = 7 



E(^Xi, Xj, X]^j^, x^ 



E{xi,xwL, 



X" 



+- 



1 

1 - ^ 



2 - Xi \ logxi 



XWh 



XWf. 



Tlogixwn) 



Xw,. 



XWh 



Ax 



Wh 



(l-f) 



2 1 {Xi Xj^{\ Xi) 

H (l - ^) 



log{xi) - 



X- 



(i-t) 



(A. 12) 



-log{xj) 



-log{x^) 
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{x^ - Xi){x^ - Xj){l - x^) 

{x.(l-f)+?(f-l) 

(1 - Xi){x^ - Xi) 

H (l - ^) 



{Xj Xi){l Xj)[X(j, Xj) 

(A.13) 



(1 - Xi)^{x^ - Xj)^ 



{x^ - Xi)^{l - x^) 



log{x^) 



B Input parameters 



Gf = 1.16 X 10~^GcV"2 , /s = 0.21 , = 5.3GeV, 
mw^ = 80.4GeV , mz^ = 91.2GeV 



log{xi) 



(A. 14) 
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